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We study numerically the propagation of two-color light pulses through a metamaterial doped with active 
atoms such that the carrier frequencies of the pulses are in resonance with two atomic transitions in the A 
configuration and that one color propagates in the regime of positive refraction and the other in the regime of 
negative refraction. In such a metamaterial, one resonant color of light propagates with positive and the other 
with negative group velocity. We investigate nonlinear interaction of these forward- and backward-propagating 
waves, and find self-trapped waves, counter-propagating radiation waves, and hot spots of medium excitation. 
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1. INTRODUCTION 

Recent advances in materials technology have led to the 
emergence of a new branch of science, the physics of 
metamaterials. These are artificially engineered compos- 
* | alexkor@math.unm.edu | XhdX exhibit properties unattainable in nature and 
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furnish exciting potential new tools for future technolo- 
giesi Electrodynamics of metamaterials has opened av- 
enues for manipulating light in unforeseen ways, as well 
as presented new paradigms in basic science. Metama- 
terials exhibiting negative refractive index within certain 
frequency ranges^^^ provide a particularly promising 
example, as they give rise to a number of unusual opti- 
cal properties impossible to observe in natural materials. 
In addition to negative refraction, i.e., reverse Snell's 
law,^^ these properties include perfect lensing^ii^ re- 
verse Doppler shift ^ reverse Cherenkov radiation / ^i-*^^ 
and, in particular, phase and group velocities point- 
ing in the opposite directions Negative refractive in- 
dex was first studied as a purely hypothetical concept 
and analyzed using simple mathematical models J^i^i^ 
Only relatively recently was this concept realized experi- 
mentally using metamaterials, first for microwaves,^ and 
subsequently for infrared and visible light. 

Nonlinearity in the interaction of light with metama- 
terials can arise both from the host material and/or 
from the embedded structures^^i^— Particularly in- 
triguing situations arise on the boundaries of negative- 
and positive-refractive-index regimes in metamaterials, 
be it in the physical or the frequency domain. The for- 
mer include, for example, nonlinear optical vortices^ 
and bistability.^^ For the latter, let us recall that nega- 
tive refractive index exists only inside limited frequency 
ranges, while outside these ranges the refractive index is 
positive. Nonlinearity can facilitate energy interchange 
among different frequency components of optical pulses 
traveling through a nonlinear metamaterial. Two possi- 
ble scenarios can occur: pulses with broad spectra whose 
wings belong to regions of different refractive-index sign, 
or multi-wave interaction in which the individual wave 
frequencies belong to such regions i^^^^^ One way 
to achieve this latter scenario is to dope a metamate- 
rial with active atoms, and launch into this metamate- 
rial two-color light pulses that resonate with two atomic 
transitions from a pair of energetically-lower level to a 
common upper level, with one transition frequency ly- 
ing in the regime of positive and the other in the regime 
of negative refraction. This transition configuration is 
known as the A configuration.-^S, Likewise, in a chiral 
metamaterial^i^ the opposite circular polarizations of 
a monochromatic light pulse that travel in the opposite 
regimes of the refractive-index sign can be coupled via 
the interaction with active atoms through a pair of tran- 
sitions in the A configuration. 

In this paper, we study two-color light pulses propa- 
gating through metamaterials containing embedded ac- 
tive atoms with two electron transitions in the A config- 
uration, corresponding to the two colors, for which one 
transition frequency lies in the positive-refraction and 
the other in the negative-refraction regime. This physi- 
cal setup provides a perfect example of nonlinear inter- 
action between wave pairs propagating within the two 
refraction regimes in the metamaterial. Since the group 
velocities of such wave pairs have the opposite signs. 



it is not surprising to find that resonantly-interacting 
waves typically travel in the opposite directions. This 
unusual interaction gives rise to a number of effects. 
The first effect we find is self-trapping, resulting in a 
pair of co-propagating solitary-wave-type pulses simi- 
lar to the self-trapped pulses in metamaterials exhibit- 
ing second-harmonic generation^^i^ The second effect 
is the counter-propagation of the linear waves shed by 
the incident pulse, traveling at the two respective car- 
rier frequencies. The third effect is the stopping of a 
pulse as it switches energy from one color to the other, 
accompanied by a hot spot of medium excitation left in 
the metamaterial at the location of the color switching 
and direction reversal. 

Several known mechanisms currently exist that pro- 
duce a negative refractive index, such as chirality,— 
hyperbolic materials and metamaterials based on 
the magnetic response utilizing plasmonic resonance ^^i^ 
All these approaches yield the same ultimate effect, and 
so, with no loss of generality, we here focus on the last 
one. We begin with Maxwell's equations for the elec- 
tric and magnetic fields, coupled to material equations 
describing the magnetic and dielectric response of the 
nanostructures and the additional response of the active 
atoms in the A configuration. We restrict our study to 
coherent phenomena, i.e., the case in which pulses are 
very short compared to all the relaxation time scales. 
Under the assumption that light pulses are short, but 
still much longer than the individual wavelengths of the 
light they contain, we have derived a set of Maxwell- 
Bloch equations that describe the electric-field envelopes 
and the medium polarizations and level-occupation den- 
sities. While we only consider idealized, lossless propa- 
gation, our equations can be readily adapted to include 
terms describing losses and pumping. This is important 
because metamaterials are quite lossy, especially in the 
negative-refraction frequency range, and doping them 
with active atoms has in fact been used to successfully 
reduce these losses^^ 

The remainder of the paper is organized as follows. 
In Sec. 12. Al we present our modeling assumptions and 
state the basic physical model that we adopt for our in- 
vestigation. In addition, we determine the fundamental 
temporal and spatial scales in this model. In Sec. 12. B| 
we sketch a derivation from this model of the Maxwell- 
Bloch equations for the envelopes of the electric-field 
components and material variables, which we use in our 
numerical simulations. In Sec. 13. A| we discuss the phys- 
ical setup underlying our numerical simulations. In par- 
ticular, we discuss the choice of the parameters, as well 
as the initial and boundary conditions. In Sees. 13. Bl and 
13. C( we discuss the phenomena we have uncovered in 
our model: co-propagating pairs of self-trapped pulses, 
counter-propagation of radiation components, and hot- 
spot formation in the medium. Finally, in Sec. HI we 
present our conclusions. 
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2. MODEL 

To set the stage for our numerical investigation, in this 
section, we briefly outline a derivation of the envelope 
equations that describe the propagation of optical pulses 
in a A-configuration metamaterial with embedded nano- 
inclusions. Ultimately, we will assume our model to de- 
scribe a metamaterial fabricated in such a way that one 
of the two allowed atomic-transition frequencies takes 
place in a regime of negative refractive index and the 
other transition frequency in a regime of positive refrac- 
tive index. Our goal is to investigate the interaction 
among waves propagating in the two regimes, and the 
resulting dynamical phenomena emerging from this in- 
teraction. 

2. A. Fundamental Equations and Metamaterial 
Model 

Electromagnetic field propagation through a dielectric 
(meta) material in the absence of macroscopic currents 
and free charges is described by Maxwell's Equations 



V X E = 



IdB 



V X H = 



V-D = 0, V-B = 0, 



(la) 



(lb) 



where E represents the electric field, B the magnetic 
induction, H the magnetic field, and D the electric dis- 
placement field. The constant c is the speed of light in 
vacuum. The interaction between the electromagnetic 
field and the (meta) material doped with active atoms is 
expressed through the constitutive relations 



B = H + 47rM, 



D = E + 47rP + 47rQ, 



(2a) 



for the magnetization M and polarization P due to the 
nanoinclusions and the polarization Q due to the active 
atoms. 

The dynamics of the magnetization M and polariza- 
tion P are modeled by the oscillator equations^"— 



47r 



E, 



"47r df^ ' 



(2b) 



(2c) 



Here, cjp is the plasmonic oscillation frequency and ujm 
is the magnetic resonance frequency associated with the 
inclusions, and 7 and j3 are form factors specific for 
the given inclusion type and material properties of the 
surrounding dielectric in the metamaterial, which can 
be evaluated by careful homogenization^^^*^ or else ob- 
tained from the measured plasmonic and magnetic reso- 
nance curves. Equation (|2bl) is the classical Lorentz 
model^ and (|2c|) is its analog in a magnetic mate- 
rial^.i^ 



The polarization contribution Q due to the active 
atoms is modeled as 



rii = 0, (2d) 



where a^, / = 1, 2, 3, are the quantum amplitudes of the 
three active atomic levels in each individual atom, J\f is 
the density of the active atoms in the material, vij are 
the electric dipole moments of the transitions between 
the l-th and j-th state, and ((•)) denotes averaging over 
all the atoms in an infinitesimally small volume at the 
location z. For every atom, the amplitudes am satisfy 
the equation 



in— -ay, 
at 



EE 

n=l 



(2e) 



for m = 1, 2, 3, which can be derived from Schroedinger's 
equation4^ Here h is Planck's constant and hwm are 
the energies of the three working levels in the active 
atoms embedded in the metamaterial, with cj^n being 
the corresponding quantum frequencies. In what is to 
follow, we will assume that they are narrowly distributed 
around three fixed values. 

For active atoms that have a pair of transitions in 
the A configuration, we assume that the frequency uji 
corresponds to the excited energy level, and U02 and uj^, 
to the two lower energy levels. In this configuration, the 
transition between the two lower levels is forbidden, and 
so the corresponding electric dipole moments vanish. 



r23 = r32 = 0. 



(3) 



2.B. Maxwell-Bloch Envelope Equations 

We assume wave propagation in the medium to be one- 
dimensional, varying in the 2:-direction only, so that the 
fields E, B, P, M, D, and H all lie in the x?/-plane. 
Equations (|lbp are then satisfied automatically. We 
non-dimensionalize equations ([T]), ()2b|) . (f2c|) . (|2d[) . and 
([2ej) using a number of fundamental scales in the sys- 
tem. We let uoq = maxj^/e \uOj — (jJk\ denote the typical 
size of the quantum frequencies uOn in Eq. ([2e|) . which 
is inversely proportional to the wavelengths of the light 
absorbed and emitted in the corresponding atomic tran- 
sitions. We let d = max |r^n| denote the typical electric 
dipole strength of the electron transitions. We define 
the cooperative frequency 



h 



(4) 



which is the typical response scale of the active atoms in 
the medium, while its reciprocal is proportional to the 
temporal width of the light pulses involved in the inter- 
action with the metamaterial. Using the scales cjq, 00 c-, 
and we rescale the spatial and temporal variables as 
(jjcz/c Uct t, the frequencies as uj/uq cj. 
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Qp = ujp/ujQ^ Qm = ^m/^o, = for 

m = 1,2,3, the electromagnetic-field and material vari- 
ables as dE/hcUc E, dR/huJc H, AirdP/huJc P, 
ATrdWL/hwc M, and the electric dipole moments as 

(Xlj =Ylj/d. 

We restrict our study to the case of the slowly- varying 
envelope approximation when light pulses include many 
electric-field oscillations, i.e., uOc <C cjq. We thus assume 
the electric field to be a sum of slowly- varying modulated 
plane waves 

E(z,t) =E+(z,t)e^^+ +E_(z,t)e^^- +c.c., (5) 

where c.c. stands for complex conjugate, uo± are the car- 
rier frequencies of the electric fields, and the fast phases 
are described by the expressions 



(6) 



We assume the fields H, P, and M to behave in a like 
manner. Here, the dimensionless wave numbers and fre- 
quencies satisfy the usual dispersion relation 



The polarization components p± , corresponding to the 
two atomic transitions in the A-configuration medium, 
are expressed in terms of the atomic level amplitudes ak 
as t^v) = i (ai(z, t)a2{z^ t)) e~*^+ and p-{z^ t, v) = 
i {ai{z,t)a^{z^t)) e~'^^- ^ where (•) denotes ensemble av- 
eraging over all the atoms in the sample in an infinites- 
imal volume at the location z and time t, whose transi- 
tion frequencies are detuned by {oOc/^o)^ from the res- 
onance with the frequencies of the respective electric 
field components. Likewise, we express the polarization 
component due to the forbidden transition between the 
two lower levels 8isp{z,t^iy) = {a2{z^t)a%^{z^t)) e^^^+~^-\ 
and the level-occupation numbers as A/" = (|ai p), n+ = 
(|a2p),and n_ = (|a3p). 

For the electric-field components Ej^ = 2E+ • a+ and 
E- = 2E_ • a_, and the macroscopic polarization and 
occupation variables, p±, p, TV, and n±, we derive the 
set of Maxwell-Bloch equations 



1 dE± 




(9a) 



k'^iuj) = cj^/i(cj)e(co'). 



(7a) 



where e(cj) and piuo) are the dielectric permittivity and 
magnetic permeability, respectively, given by the expres- 
sions 



e{u) = 1 + 



7^1 



M ■ 



which follow from the oscillator equations (|2b[) describ- 
ing the medium polarization P and magnetization M. 

For the polarization Q due to the active atoms, de- 
fined in Eq. (|2dp . we assume the form 



/oo 
[p+{z,t,iy)ocle'''+ 
-OO 

i/)ale^^- - c.c] g{u)du, (8) 

where q:+ = cxi2 and a_ = a 13, and the components 
p±{z,t,u) are detuned from the exact resonance with 
the carrier frequencies uj± of the electric-field compo- 
nents E±(z,t) by the amount (cJc/^o)^- In other words, 
we assume that the respective atomic transition fre- 
quencies Qi — Q2 and — Qs are related to and 
uj- by the relations cj^ = l^i — 1^2 + (uJc/cjOq)!^ and 
uj- = 1^1 — 1^3 + (cjOc/^o)^' The detuning u arises mainly 
due to the motion of the active atoms, which gives rise 
to a Doppler shift in the transition frequencies; u is 
therefore distributed according to a function g{u) with 
J^oo 9{^)dh' = 1, which describes the shape of the spec- 
tral line.^^ In assuming that the actual detuning is of 
size (cJc/cJo)^, we assume that this line is very narrow as 
compared to the typical size of the variations of the di- 
electric permittivity e(a;), magnetic permeability 
and the wave vector k{uj) with the frequency uo. 



^=iyp^+'^[E+{N-n+)-E.p*l (9b) 



dp- 
dt 



ivp- 



■n.)-E+p], (9c) 



■E-pX]. 



(9d) 



dN 
'dt 



-[E+pl + EXp+ + E_p*_ + E*_p_], (9e) 



'-^ = \[E,pl^Elp^], 



(9f) 



where p± = p{oJ±) and e± = e{uj±)^ and the group ve- 
locities v± = v{uj±) are computed from the dispersion 
relation ([7a|) via the formula v{uj) = l/k\uj)^ with the 
prime denoting the derivative with respect to the fre- 
quency cj, which leads to the expression 




1^% 



(10) 



M 



We will investigate the solutions of equations ([9]) numer- 
ically in the remainder of this paper. 
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3. RESULTS 

Our goal in the remainder of this paper is to investigate 
a doubly-resonant interaction between light pulses and 
a metamaterial doped with active atoms in the A config- 
uration, such that one of the two corresponding atomic 
transitions takes place in the positive and the other in 
the negative refractive-index regime. Using numerical 
simulations, we uncover several phenomena that arise as 
consequences of this interaction. 

3. A. Physical Setup 

We simulate light pulses with carrier frequencies which 
are centered around the values uoj^ and cj_, as depicted 
in Fig. [11 around which the transition frequencies of the 
atoms in the A configuration embedded in the resonant 
metamaterial, l^i — Q.2 and l^i — 1^3, are also narrowly 
centered. In a typical metamaterial, light interacts reso- 
nantly with metallic nanoinclusions embedded in it, and 
due to the resonant nature of the negative refraction, 
the value of the refractive index is negative, n < 0, only 
within a limited interval of frequencies of uj. We will 
study the case when cj+ is in the regime of positive re- 
fraction, n > 0, and uo- is in the regime of negative 
refraction, n < 0, see Figure O so that the group ve- 
locities of the two pulse components point in opposite 
directions. To best focus on the opposite directional- 
ity of the two electric-field-component group velocities 
v±^ we consider an idealized version of equations (|9]) in 
which we set all constants to 1 , except v± , which we set 
to v± = ^1. This choice may also correspond exactly to 
the case of a chiral metamaterial, in which the two A- 
configuration transitions in the embedded active atoms 
correspond to the same frequency but different polar- 
izations of light J ^^1^^ In addition, we consider limit of a 
narrow spectral line, when g{v) = S{i'), where S{-) is the 
Dirac delta function. 

The physical setup that we consider in our simula- 
tions is as follows: Initially, the system is in the ground 
state, N{z, 0, = 0, which is characterized by a relative 
distribution of the ground- level populations n+(z,0,i^) 
and n_ (z, 0, z/); the medium polarization due to the one- 
photon transitions is absent, p+(z,0,z^) = p_(z,0,i^) = 
0; the ground states are coupled through the initial po- 
larization p{z^0^iy) = po{z^iy). There is no light in 
the medium initially, £^+(z,0) = E-{z^O) = 0, and a 
right-propagating incident pulse, £^+(0,t) = £^+o(^), is 
injected at the left edge of the metamaterial sample. 
We set E-{l^t) = 0, where / is the 2:-coordinate of the 
right edge of the sample, as we inject no left-propagating 
electric-field component into the metamaterial. In the 
simulations to follow, we take for the boundary condi- 
tion Ej^oif) the Gaussian function given by 

E+oW = Ae"^'"''^^'/'"', (11) 

where A is the amplitude, td is the delay time, and a the 
width of the injected pulse. Specifically, for the latter 
two parameters, we fix the dimensionless quantities td = 



3 and a = 1 throughout our investigation. In addition, 
we fix the length of the sample to be / = 40. 

In the majority of previous studies, the initial polar- 
ization component po(^, ^) was assumed to vanish, which 
corresponds to the absence of coherence between the two 
lower-energy states in the A-configuration medium. A A- 
configuration medium with nonzero po can be prepared 
by irradiating it with a two-color electromagnetic field in 
such a way that the transition between the two lower lev- 
els is in a two-photon resonance. This can be achieved if 
the difference between the frequencies of the two colors of 
the light used in preparing the medium equals cj+ — a;_ . 
Since a dipole transition between these two lower lev- 
els is forbidden, the induced coherence between them, 
refiected in the polarization component po, will persist 
much longer than other polarization components. The 
nonzero component po provides an immediate, nonlinear 
coupling between the two electric field components E±^ 
regardless of their size. 




Fig. 1. A-configuration system with transition frequencies cj^ 
and UJ- . The transition with the frequency 0;+ corresponds to 
the response of the metamaterial in a regime with a positive 
index of refraction, and the transition with the frequency U- 
corresponds to the response of the metamaterial in a regime 
with a negative index of refraction. The medium is initially 
prepared so that a nonzero component of the polarization, 
Po, corresponding to the two-photon transition between the 
two lower levels, is present. 

We are primarily interested in the physical setup in 
which a right-traveling electric field component is in- 
jected into a metamaterial containing active atoms in the 
A configuration, and a left-traveling electric field compo- 
nent is generated via color switching. We are also inter- 
ested in exploring the precise details of how this second 
color of light propagates in the opposite direction of the 
original input light. A schematic of this setup situation 
is displayed in Figure [H Here, the red light, labeled 
pump, is injected into the metamaterial and propagates 
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n>0 n<0 



|E^(x,t)| 



CO. 



'+ 



co_ 



Fig. 2. Due to the resonant nature of the negative refractive 
index phenomenon, the sign of the refractive index depends 
on the frequency. The intermediate (forbidden) range of fre- 
quencies, in which e and /a have (real parts with) opposite 
signs and there is no wave propagation, is depicted by the 
green dot. 



Pump 




A— Medium 



Fig. 3. A schematic of the physical setup investigated numer- 
ically: in this case only the red light, the pump is sent into 
the A- configuration metamaterial. This red light is partly 
converted into blue light, which travels in the opposite direc- 
tion. 



to the right. The goal is for the red light to be converted 
into blue light, labeled signal, which would propagate in 
the negative refractive index regime and so to the left, 
and eventually exit the medium at the same end at which 
the red light entered it. 

We now describe the phenomena we have observed in 
our numerical simulations, which arise as a consequence 
of the color conversion discussed above. 

3.B. Self- Trapping and Counter-Propagation 

The results of the simulations we display in Figs. H] and 
[5] show a self-trapped, co-propagating pair of nonlinear 
pulses with oscillating amplitudes, which forms in both 
electric-field components, and from the incident 
pulse with the nonzero component E_^Q{t) in Eq. ([TT]) . 
These components travel simultaneously in the same di- 
rection and with the same velocity despite the fact that 
the linear group velocities v± corresponding to their re- 
spective refraction regimes point in the opposite direc- 
tions. 

For the simulation depicted in Fig. |4] , the amplitudes 
E-^ and E- of the self-trapped pulse pair contain virtu- 
ally all the energy of the incident pulse. These ampli- 
tudes are almost equal, and the resulting trapped pair of 
pulses appears to propagate at a constant velocity, and 
considerably more slowly than the group velocity v+ of 
the positive refraction regime. The two co-propagating 
pulses maintain rather complicated, in-phase, oscilla- 
tory shapes during their propagation, and shed a small 
amount of radiation. 

In Fig. O we display a different setup in which the 
amplitude of the incident pulse E-^o{t) and the initial 
polarization po are smaller, while the initial population 




5 10 15 20 25 30 35 40 
|E.(x,t)| 




5 10 15 20 25 30 35 40 

m 

Fig. 4. Self-trapped, co-propagating pulse pair. Non- 
vanishing initial conditions: po = 0.3, n+ = 0.8, n- = 0.2. 
Amplitude of the injected Gaussian pulse E-^o{t): A = 2. 



difference has the opposite sign from those used to gen- 
erate Fig. m In this case, the emerging co-propagating 
pulse-pair not only propagates more slowly than the lin- 
ear group velocity v+ of the positive refraction regime, 
but also appears to decelerate during propagation. Note 
that the amplitude of the pulse-component E- shown in 
in Fig. [5l propagating in the negative refraction regime, 
is considerably smaller than that of the pulse-component 
propagating in the positive refraction regime. 

The spatio-temporal distribution of the electric-field 
component in Fig. [5] shows a forerunner wave, prop- 
agating with the group velocity V-^ positive refraction 
regime and with an oscillatory distribution of its tails, 
followed by the self- trapped pulse, and the waves in the 
region behind this pulse. The forerunner wave creates 
an oscillatory continuous-radiation field in the compo- 
nent E- , with a much smaller amplitude than the com- 
ponent E-^ , which both propagate in the region between 
the forerunner wave and the self- trapped pulse. Another 
source of continuous radiation is the self-trapped pulse 
itself, which largely sheds radiation in the region behind 
it. Note that these two regions of continuous radiation 
exhibit very distinct features: One can easily notice the 
radiation in E- emerge from the peaks in the oscillations 
of the self-trapped pulse and propagate backward in the 
region behind this pulse, as indicated by the arrow in 
Fig. [5l Note that this radiation still travels more slowly 
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than the group velocity V- in the negative refraction 
regime. 

From the middle two panels in Fig. O one can clearly 
see the excitation of the medium polarizations p+ and 
p_, induced by the co-propagating pulse pair. A closer 
look shows that the electric-field and polarization oscil- 
lations along this pulse have the opposite phase, indicat- 
ing that an oscillatory exchange of energy between the 
electric field and the medium takes place as the pulse 
propagates. Note the sharp dip of the polarization com- 
ponent p+ along the pulse trajectory, which is where 
the amplitude of this component changes sign. In addi- 
tion, a depletion of the electron population density n+ 
appears along the pulse trajectory, accompanied by a 
corresponding increase in the population density n_. 

3.C. Hot Spot 

The incident pulse with the nonzero component Ej^ , dis- 
played in Fig. [6] propagating in the positive-refraction 
regime, slows down to a halt and decays in time until it 
is virtually extinguished. Part of this pulse switches into 
the component E- , which contains a great deal of radi- 
ation propagating backward. Note that the amplitude 
of the incident pulse Ej^^it) is the same as in Fig.|H the 
initial polarization po is the same as in Fig. [5l and the 
initial population densities are equal. 

In addition, we see that at the location at which the 
entire incident pulse has disintegrated and switched into 
the component a strong excitation in the medium 
polarization p+ and depletion of the occupation density 
n+ emerges. A considerable amount of energy appears 
to be deposited in the medium at this location. This 
energy forms a stationary hot spot of medium excita- 
tion, which also does not appear to noticeably decay in 
time. The polarization p_ and occupation density n_ 
appear to also refiect the underlying backward propa- 



gating continuous radiation in the E- component of the 
electric field. 

4. CONCLUSIONS 

In this paper we have investigated the nonlinear cou- 
pling and consequent flow of energy between light pulses 
propagating in the positive refractive index regime and 
light pulses propagating in the negative refractive index 
regime (or vice versa) within a resonant metamaterial 
doped with active atoms in the A configuration. 

The phenomena we have observed in our numerical 
simulations include the creation of a self-trapped, co- 
propagating, nonlinear pulse pair that propagates with 
neither of the group velocities. In particular, this pulse 
behaves in an oscillatory fashion and propagates in the 
same direction but with a velocity slower than the orig- 
inal incident pulse. The radiation this pair sheds typi- 
cally propagates in the directions indicated by the linear 
group velocities of its respective color. 

We have also observed a pulse with the color corre- 
sponding to a positive-index regime to almost completely 
switch into the negative index regime where radiation is 
created, which travels backward. Simultaneously, a hot 
spot of permanent medium excitation is formed when 
the pulse in the positive index regime slows and eventu- 
ally stops. From the location of this hot spot, the electric 
field in the negative index regime sheds backward prop- 
agating radiation. The energy deposited in the hot spot 
appears to remain stationary in space as time increases. 

In the simplest, most idealized case of light interact- 
ing with a A-configuration optical medium and in the 
absence of any initial medium polarization excitation, 
we have predicted a particularly transparent version of 
polarization and/or color switching In particular, if 
one of the two energetically-lower levels is initially more 
populated along the medium sample than the other, the 
electric-field polarization or color of any pulse propagat- 
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ing through this medium switches until the pulse inter- 
acts with the medium only through the transition be- 
tween the excited upper level and the less populated of 
the two lower levels. In the presence of a non- vanishing 
initial polarization po i the system instead asymptotically 
approaches a specific combination of the two transitions, 
whose details depend on the value of Further sim- 
ulations aimed at drawing analogies and comparisons 
between these idealized results and the dynamics of the 
light-metamaterial interaction investigated here will be 
conducted in future investigations. 
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